Cold heteronuclear atom-ion collisions 
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We study cold heteronuclear atom ion collisions by immersing a trapped single ion into an ultra- 
cold atomic cloud. Using ultracold atoms as reaction targets, our measurement is sensitive to elastic 
collisions with extremely small energy transfer. The observed energy-dependent elastic atom-ion 
scattering rate deviates significantly from the prediction of Langevin but is in full agreement with 
the quantum mechanical cross section. Additionally, we characterize inelastic collisions leading to 
chemical reactions at the single particle level and measure the energy-dependent reaction rate con- 
stants. The reaction products are identified by in-trap mass spectrometry, revealing the branching 
ratio between radiative and non-radiative charge exchange processes. 

PACS numbers: 34.50.Cx,34.70.-|-e, 37.10.Ty 



Cold collisions are characterized by the de-Broglie 
wavelength of the colliding particles becoming compa- 
rable to the length scale of the molecular interactions. 
Quantum mechanics then dominates the elastic and in- 
elastic scattering phenomena and cross sections. Under- 
standing elastic collisions in this regime is fundamen- 
tal for achieving and optimizing sympathetic cooling of 
species which do not allow for direct laser cooling, for ex- 
ample complex molecules. On the other hand, inelastic 
collisions of translationally cold reaction partners provide 
a unique opportunity to observe quantum mechanical de- 
tails of chemical processes which typically are averaged 
out at higher temperatures p|. It is expected that the 
detailed understanding of cold collisions will eventually 
pave the way towards coherent chemistry, in which full 
quantum control of the internal states and reaction path- 
ways will be possible. 

Over the past few years, impressive advances have 
been made in exploring cold collisions of neutral atoms. 
This has lead, for example, to the controlled creation of 
weakly bound dimers [4I and trimers 3| and the produc- 
tion of molecular gases in the rovibrational ground state 
[^. In comparison, cold ion- neutral collisions are much 
less explored. Ion-neutral reactions offer the outstand- 
ing experimental possibility that reactants and reaction 
products can be observed and manipulated at the sin- 
gle particle level and can be trapped for very long times 
allowing for precise measurements 'sM!!!. Collisions of 
cold (millikclvin) trapped ions with hot (~ 10^ K) neu- 
tral atoms have provided insights into chemical reactions 
[a,|7|,|lO| and play an important role in understanding the 
molecular composition of the interstellar gas |11]. Only 
very recently, the regime in which both collision part- 
ners are translationally cold has been accessed by us- 
ing trapped ions interacting with cold molecular beams 
[7| or trapped neutral atom samples 12|, |l3[ . Neverthe- 
less, elastic scattering in the quantum regime of many 
partial waves with its decisive energy dependence of the 
cross section [a oc E^^^^) was not yet observed. Regard- 
ing charge exchange, experimental observations 12] were 



well explained by the Langevin model (ctl oc i?^^/^) that 
is based on classical mechanics to describe the mobility 
of ions in gases (l4| . 

In this manuscript, we explore heteronuclear collision 
processes between trapped atoms and ions in which both 
collision partners are translationally cold. We study elas- 
tic and reactive collisions using a single ion only in or- 
der to better control the kinetic energy than in larger 
ion crystals. This enables us to investigate the energy 
dependence of elastic and inelastic scattering processes 
from 20 to 450 /^eV (equivalent to thermal energies in 
the range of 0.2 to 5 K). We find the elastic collision rate 
significantly larger than predicted by Langevin theory 
but in full agreement with a quantum mechanical calcu- 
lation. Moreover, we investigate inelastic collisions and 
measure the energy dependence of the charge exchange 
cross section and the branching ratio for radiative vs. 
non-radiative charge exchange. This provides a compre- 
hensive characterization of binary atom-ion collisions in 
this energy range. 

Atom-ion scattering is dominated by the polarization 
potential, which asymptotically behaves as V{r) = — |^, 
and a hard-core repulsion at the length scale of the Bohr 
radius. Here r is the internuclear separation, C4 — 
q denotes the charge of the ion, a is the dc polar- 



izability of the neutral particle, and eg is the vacuum per- 
mittivity. Atom-ion scattering was theoretically investi- 
gated as early as 1905 when Langevin calculated the drift 
velocity of ions in a buffer gas. The Langevin model is 
based on classical mechanics and the predicted cross sec- 
tion is ctl = TT^y2Ci/E, leading to an energy independent 
collision rate constant. Langevin-type collisions occur in 
close encounters between atom and ion when the impact 
parameter is below a critical value be = {2C4/Ey^'^ [l9| 
and exhibit a large momentum transfer between atom 
and ion. Going beyond the Langevin model, the quan- 
tum mechanical description of scattering yields the an- 
gular dependence of the differential cross section in more 
detail. In the energy range under consideration, elastic 
scattering from the polarization potential involves sev- 



eral partial waves since the s-wave scattering regime of 
atom-ion collisions is only at sub-microkelvin tempera- 
tures [iSl |l6[. The semi-classical approximation of the 
full quantum mechanical elastic scattering cross section 

for a collision energy E is a{E) = 7r(l + 7r^/16) ( j^ j 

la lul (m denotes the reduced mass). This includes in 
particular a large contribution towards forward scatter- 
ing under very small angles [l8|, mainly from the cen- 
trifugal barrier. The forward scattering peak results from 
interference of partial waves similar to the bright 'Pois- 
son spot' emerging behind a beam block in wave optics. 
In these collisions the momentum transfer between the 
collision partners is extremely small and very difficult to 
detect. Until now, neither cold ion-neutral collisions nor 
ion-mobility spectrometry could so far access this regime. 
In our experiment we study collisions between ultra- 
cold ^^Rb atoms and a single Yb+ ion. We prepare 
2.2 X 10® neutral atoms in the \F — 2, mp — 2) hyper- 
fine ground state at T ~ 300 nK in a harmonic mag- 



netic trap of characteristic frequencies (uj. 



27r X (8, 28, 28) Hz [13|,|20|. We spatially overlap the neu- 
tral atoms with a single Yb"*" ion trapped and laser cooled 
in an RF-Paul trap. The ion in the Paul trap expe- 
riences a pseudo-potential which is derived from time- 
averaging a very rapidly oscillating electric quadrupolc 
field E(a;, y) = ^^^T (x, —y, 0). Vrf denotes the ap- 
plied voltage, R is the distance of the ion to the elec- 
trodes, and J7 = 27r X 42.7 MHz is the drive frequency. 
The time-averaged potential provides a harmonic con- 
finement for the ion with characteristic trap frequencies 
wj_ = 27r X 150 kHz. On top of the time-averaged motion, 
the ion undergoes very rapid oscillations in position (mi- 
cromotion) of an amplitude amm which is proportional 
to the local electric field strength. The micromotion is a 
driven oscillation in which the ion's energy is on average 
conserved because momentum of the ion and the applied 
force are 7r/2 out of phase. Axial confinement of the ion 
with a characteristic frequency of ujax = 27r x 42 kHz is 
provided by electrostatic fields. 
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FIG. 1. (Color online) a) Time-resolved fluorescence of a 
single ^^^Yb ion. The signal is averaged over 69 repetitions 
and the error bars are determined by the photon shot noise. 
The line is a fit [2l|] resulting in a mean energy of the ion of 
(230 ± 9) fieV. b) Measured mean ion energy vs. offset elec- 
tric field used for tuning the excess micromotion. The line is 
a linear fit to the data. 



A collision with an ultracold atom changes the instan- 
taneous velocity of the ion and disrupts its coherent mi- 
cromotion oscillation. As result, energy can be trans- 
ferred from the ion to the neutral atoms, leading to cool- 
ing of the secular motion of the ion [13[ , but also energy of 
up to ^ ^21^ ^2^2^^ ^g^^ |-|g transferred from the driving 



field to the ion's secular motion 22|-j24] . Depending on 
the exact parameters, the ion will equilibrate at a mean 
energy determined by the two processes. We tune the 
energy of the atom-ion collisions by applying a homoge- 
neous offset electric field Eq transverse to the symmetry 
axis of the ion trap. This displaces the ion from the 
geometric center of the Paul trap and introduces excess 
micromotion of amplitude amm.ex = \/2qEo/{mion<jJ±^)- 
rriion is the mass of the ion. We center the neutral atomic 
cloud at the new equilibrium position by applying homo- 
geneous magnetic fields of ^ 100 mG. 

We investigate the effect of the excess micromotion 
o.mm,ex On the equilibrated mean energy of the ion after 
it has interacted with a neutral atom cloud for <o = 8 s. 
During the interaction with the neutral atoms the ion 
is in its electronic ground state. For the energy mea- 
surement we release the neutral atoms from the trap and 
set the offset electric field to zero. Then, we illuminate 
the ion with laser light red-detuned by A = — 0.25F 
from the S'1/2 ^ -P1/2 transition of Yb+ at 370 nm, 
(F = 27r X 20 MHz denotes the linewidth of the atomic 
transition) and monitor the temporal increase of the flu- 
orescence rate (see Fig. [1^). We fit the data with the 
solution of the time-dependent optical Bloch equations 
for the motion of the particle in the pseudopotential and 
perform a thermal ensemble average 2l|,|25J. In Fig. [TJd 
we display the measurement of the mean energy of the 
ion vs. the applied offset electric field Eg which exhibits 
the expected quadratic dependence. We use this energy 
tuning mechanism to study the energy dependence of the 
collision cross sections in the next paragraphs. The min- 
imum energy is determined by residual micromotion. We 
find that the ion equilibrates after less than 100 ms inter- 
action time to its average kinetic energy. 

The kinematics of a binary collision between an ion 
of energy Eion and a neutral atom of mass m„, which 
we assume at rest (£'„ — 0), results in a relative en- 
ergy transfer l^EjE^on = 47/(7 -h 1)^ sin2(6'/2) with 
7 — rriion /"mn- Collisions in which the scattering angle 9 
exceeds 9 cut = 2aicsvi\{^Ecut{l + lY /{-i^Eion)) lead to 
atom loss from the trap. Ecut is the effective trap depth 
of the atom trap. A neutral atom after a collision will 
be detected as lost in absorption imaging if it moves at 
a very large orbit in the trap and does not thermalize 
with the neutral atom cloud within the time ig- Numeri- 
cal calculations show that thermalization can be achieved 
only for energy transfers less than 0.8 neV. In Fig. [2^ we 
show the measured atom loss as a function of the ion 
energy. We find that the atom loss exceeds the classi- 
cal Langevin prediction AA'^ = riQcrLvto by more than an 



order of magnitude and, moreover, shows a pronounced 
energy dependence, also not explained by the Langevin 
model. Here Hq is the local neutral atom density and v 
is the relative velocity of atom and ion. To model our 
data, we numerically calculate the quantum mechanical 
differential elastic scattering cross section I{9,E). We 
calculate the scattering phase shifts in a semi-classical 
approximation for large values of the angular momentum 
for a —Ci/2r'^ potential and perform a random-phase ap- 
proximation for small angular momenta for the unknown 



short-range repulsion [15[. Using I{9,E) we calculate 



the expected neutral atom losses in a Monte Carlo simu- 
lation. The solid line in Fig. [2^ shows the predicted loss 
for large energy transfer {/S.E > Ecut)- 
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FIG. 2. (Color online) a) Atom loss from the magnetic trap vs. 
mean energy of a single ^'''^Yb^ ion. The experimental data 
are averaged over ~ 70 realizations each and the standard 
error is given. The solid line is a numerical simulation based 
on the full elastic scattering cross section. The dashed line 
shows the prediction of the Langevin model, b) Temperature 
increase of the neutral atom cloud vs. mean energy of a single 
^^^Yb^ ion. The solid line shows the numerical results using 
the full elastic scattering cross section. The dashed line shows 
the maximally possible evaporative heating effect. The mean 
energy of the ion is derived from Fig. [T] 

In the case of small energy transfer [0 < 9 cut) the atom 
remains in the trap and thermalizes with the remaining 
atom cloud. The very small amount of transferred energy 
and momentum {Ocut is in the mrad-range) make these 
collisions very challenging to observe. However, using 
ultracold atoms at nanokelvin temperatures as reaction 
partners allows us to resolve a relative energy transfer 
on the order of ~ 10~^, inconceivable in previous exper- 
iments. The observed temperature increase (Fig. ^Bp) is 
explained by two different contributions. Small energy 
transfers {/S.E < Ecut) are dominant. The other effect 
results from selectively removing atoms from the center 
of the cloud where their potential energy is below the 
thermal average, leading to 'evaporative heating'. The 
dashed line shows the maximally possible contribution 
from evaporative heating, assuming all observed losses 
of Fig. [2^ being from the center of the trap. The solid 
line is the predicted temperature increase taking into ac- 
count both contributions. The simulations in both Fig. 
[5^ and Fig. ^jp use the same value of Ecut as the only free 
parameter and we obtain Ecut = 0.8 neV, in agreement 



with the experimentally expected value. 

Reactive chemical processes in atom-ion collisions, like 
charge exchange or molecule formation, can only occur 
when atom and ion approach each other to distances at 
which the electronic wave functions overlap, which is the 
case in close encounters. In the regime of many partial 
waves, the rate constant for chemical processes is thus 
predicted to be proportional to the Langevin rate con- 
stant aLV and to be independent of the collision energy 
la . Iig . l26j . For equal chemical elements strong loss due 
to near-resonant charge exchange was observed [I4I in 
agreement with the Langevin estimate. 
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FIG. 3. (Color online) a) Schematic of the molecular poten- 
tials, b) Measured ion loss probability as a function of offset 
electric field. The ion's mean energy is derived from Fig. [T] 
Each data point is averaged over approximately 100 repeti- 
tions of the experiment. The dashed lines show the average 
of the data, c) Position dependence of the ion loss probabil- 
ity for ^^^Yb"*" as the ion is scanned through the cloud. The 
red line is the theoretical density profile of the thermal cloud. 
d) The linear dependence of the ion loss rate of the density 
signals binary collisions leading to charge exchange reactions. 
Each data point is averaged over approximately 350 repeti- 
tions of the experiment. 

In our case, the Rb (5s ^S'1/2) atom in |F = 2, mp = 2) 
and the Yb"^ (Gs^S'i/a) ion in \J = 1/2, mj = ±1/2) 
collide in the excited A^E+ singlet or the a^S+ triplet 
state of the (RbYb)"*" molecular potential (see Fig. [3^). 
The electronic ground state X^T,^, which asymptotically 
corresponds to Rb+ (ip^^So) and neutral Yb(6s^^5'o), 
lies 2.08 eV below. The most striking experimental ob- 
servation of a chemical process is the disappearance of 
the Yb+ fluorescence which we detect after the interac- 
tion with the neutral atoms. Using the same ion energy 
tuning mechanism as described above, we have investi- 
gated the energy dependence of the reaction rate con- 
stant. For two different isotopes, ^^^Yb and ^^^Yb, we 
find the rate constant to be independent of the collision 
li,[ll] (see Fig. ^). The average 



(2.8 ± 0.3) X 10-20 m^/s and 



energy, as predicted 
rate constants of -R172 
_Rj74 = (4.0 ± 0.3) X 10^2" m'^/s are five orders of mag- 
nitude smaller than in the homonuclear case jl2| and of 
the same order of magnitude as predicted for the similarly 



asymmetric system (NaCa)"*" [26|. There is a systematic 
error between the measurements for the two isotopes of 
15% due to uncertainty of the density determination of 
the thermal cloud. 

Additionally, we have measured the density depen- 
dence of the inelastic atom-ion collisions by position- 
ing the ion at different locations inside the neutral atom 
cloud (Fig. [3]:). We achieve this without increasing the 
ion's micromotion by applying a homogeneous magnetic 
bias field to displace the position of the neutral atomic 
cloud. Monitoring the ion loss rate as a function of the 
position, we find that the inelastic atom-ion collisions 
scale linearly with the local atomic density, indicating 
that the charge exchange reactions are binary atom-ion 
collisions (see Fig. [3li). This also demonstrates the prin- 
cipal capability of the single ion as a local probe for den- 
sity measurements in a neutral atom cloud [27[. 
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FIG. 4. (Color online) a) Mass spectrometry signal together 
with dashed lines indicating the expected frequencies. The 
black trace corresponds to two ^^^Yb"'' ions. For one Rb^ 
ion and one Yb^ ion (red curve) the collective mode fre- 
quency is higher. No signals at the resonance for Yb"*" and 
(RbYb)"'" were obtained, b) Observed distribution of the re- 
action products for two Yb^ ions after 8 s interaction time 
with the neutral Rb atoms (blue bars, 486 events). The red 
bars (543 events) show a comparative measurement without 
neutral atoms. 

Charge exchange can occur by emission of a photon 
(radiative charge exchange) or as a nonadiabatic tran- 
sition between molecular levels. For low temperatures, 
radiative charge exchange has been predicted to be the 
dominating process for (NaCa)"*" [26|. In order to investi- 
gate the reaction products of the charge exchange process 
we perform mass spectrometry. To this end, we load two 
Yb"*" ions into the ion trap and overlap them with the 
neutral cloud. In the cases in which only one of the two 
ions undergoes a reaction, the other ion serves to identify 
the reaction product by measuring a common vibrational 
mode in the ion trap [3, 01 ■ We excite the axial secular 
motion of Yb^ by applying intensity modulated light at 
370 nm at an angle of 60° relative to the symmetry axis 
of the Paul trap and monitor the fluorescence. The fre- 
quency of the intensity modulation is swept linearly from 
30 kHz to 55 kHz in 2.5 s. When the intensity modulation 
coincides with a collective resonance of the secular mo- 
tion, the ions are heated and the fluorescence reduces 
(see Fig. 2^). The collective mode for a Rb+ and an 



Yb+ ion in the trap is 13% above the bare Yb"*" mode at 
42 kHz, whereas the mode of (RbYb)"*" and Yb'^ would 
be 12% lower in frequency. Fig. ^p shows the histogram 
for the observed charge exchange processes. In the cases 
in which one Yb"*" ion is lost, we find approximately 30% 
probability for the production of cold Rb+ and 70% for 
a complete loss. If the reaction products take up 2eV as 
kinetic energy in a non-radiative charge exchange they 
would be lost due to finite depth (~ 150 meV) of our ion 
trap. We have not observed the formation of a charged 
molecule in this process. 

In conclusion, we have provided a comprehensive sur- 
vey of cold heteronuclear atom-ion collisions. We have 
observed that Langevin theory is insufficient to describe 
cold atom-ion collisions but the full quantum mechanical 
cross section is required to describe our data. Moreover, 
we have observed charge exchange reactions between a 
single ion and ultracold atoms and analyzed their prod- 
ucts. Our results provide an excellent starting point for 
future experiments targeting the full quantum control of 
chemical reactions at the single particle level. Photoas- 
sociation or magnetic field induced Feshbach resonances 
[la | could be used to create single trapped cold molecules 
in specific rovibrational states. 
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